Mon. Not. R. Astron. Soc. 000,[T]g](2007) Printed 3 May 2008 (MN MbX style file v2.2) 



An accretion disc model for quasar optical variability 



Shuang-Liang Li*, Xinwu Caof 

Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, 200030 Shanghai, China 



Accepted 2008 April 10. Received 2008 April 3 



ABSTRACT 

Some different correlations between optical-UV variability and other quasar properties, 
such as luminosity, black hole mass and rest-frame wavelength, were discovered. The positive 
correlation between optical-UV variability amplitude and black hole mass was first found by 
Wold et al., and this was confirmed by Wilhite et al. We suggest that the accretion disk model 
can explain these correlations, provided the optical-UV variability is triggered by the change 
of accretion rate. The disk temperature of accretion discs decreases with increasing black hole 
mass, which leads to systematical spectral shape difference with black hole mass even if the 
black hole is accreting at the same rate m (to = M/ MEdd)- The observed positive correlation 
between optical-UV variability and black hole mass can be well reproduced by our model 
calculations, if the mean accretion rate toq ^ 0.1 with variation of dm ^ 0.4 — O.Stoq- We 
also found that the observed correlations of optical-UV variability amplitude with luminosity 
or rest-frame wavelength can be qualitatively explained by this accretion disc model. 
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1 INTRODUCTION 

' Quasars are variable in different wavebands, and their vari- 
, ability has been extensively explore d for sev eral decades 

■ ( ISmith&Hoffleit| [r963: Matthews & San dage|[T963h . It has been 
accepted that quasars are powered by black hole accretion, and their 

■ variability therefore provides an important clue to understand the 
' physics at work in accretion discs. The timescales of optical vari- 
ability for quasars vary from hours to decades. The variability is 
found to be correlated with other quasar properties, such as redshift 
and luminosity. There is a strong negative correlation of variability 
with quasar luminosity, while a significant p ositive correlation be- 
tween variability and red shift is also present jCristiani et al.|[l993 : 
IVanden Berk et alj|2004l) . An anti-correlation between variability 
amplitude and rest-frame wavelength has also been observed, so the 
correlation between variability and redshift was explained by the 
selection effect, because quasars at higher redshifts are more ea sier 
to be probed in shorter rest -frame wavel ength t Cid Fernandes et al.l 
1 19961 : ICristiani et alJ[T996 ). However, Vanden Berk et afj l l2004h 
claimed that the evolution of quasar population or the change of 
quasar variability mechanism can cause quasars to be more vari- 
able at higher redshifts based on a sample of over 25,000 quasars 
from the Sloan Digital Sky Survey, which is different from the se- 
lection effect s cenario. 

Recently, IWold et al.l ilOOH) matched quasars from the the 
Quasar Equatorial Survey Team Phase 1 (QUEST 1) variability sur- 
vey with broad-line objects from the Sloan Digital Sky Survey, 
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which leads to ~ 100 quasars. The black hole masses of these 
quasars are estimated from the widths of broad-line Ha. They mea- 
sured the variability amplitudes from the QUEST 1 light curves, 
and found a correlation between quasar variability amplitude and 
black hole mass. They found that the correlation seems to be un- 
likely caused by selection effects or the contamination coming from 
host galaxy or other well-known correlations. This correlation be- 
tween black hole m ass and variability amplitude was confirmed by 
IWilhiteetal.1 ( l2007l) by using a sample over 2500 quasars. 

There were many efforts devoted to reveal the physics 
behind the correlations between quasar optical-UV variability 
and other quasar properties. Many different scenarios were sug- 
gested, such as, variation of accretion rates, d isc instabilities, stel- 
lar collisions, a n d mic r olensing effect (e . g., IPerevra et alj|200^: 
'Mineshige et al.' 'l994'; 'Kawaguchi et al.' 'l998"; 'Terlevich et alj 
1992; TorriceUi-Ciamponi et al. 2000; Cid Fernandes et al. 200^ 
[Morgan et al.ll2006f) . It is widely believed that quasars are powered 
by black hole accretion. However, the spectral shape in optical-UV 
bands for a standard thin accretion disc can be approximated as 
fi, oc u^^^, while the observed quasar continua usually exhibit 
oc ly^^. Gaskell (2007) suggested that the observed quasar spec- 
tra can be reproduced by accretion discs with a temperature gradi- 
ent of r(_R) oc instead of r(i?) oc i?"" '^^ as predicted by 
standard thin disc model. This means that there is some flow of heat 
outwards in the disk. They use the optical variability timescales of 
quasars to argue that the variations must propagat e at close to th e 
speed of ligh t, rather than on viscous timescales ( lGaskellll2007h . 
IPerevra et al.l (|2006) found that the multi-epoch quasar spectra in 
the rest-frame wavelength range 1300-6000A can be fitted by the 
standard thin disc model, provided their accretion rates vary from 
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one epoch to the next. Thus, it is a reasonable scenario that the 
variability may be caused by the change of accretion rates. 

In this work, we investigate the variability amplitude of 
quasars as a function of black hole mass based on accretion disc 
models assuming the variability to be caused by the change of ac- 
cretion rates. We also compare our theoretical calculations with the 
observed correlations of variability amplitude with other quantities 
(i.e., luminosity or observed wavelength). 



2 ACCRETION DISC MODELS 

The big blue bump observed in many AGN is believed to be re- 
lated to some kind of thermal emission (Malkan & Sargent 198^; 
IShalvapin et"ai]|2002l) . The thermal emission of the standard thin 
accretion disc can reproduce the general features of AGN continua, 
provided some other components, such as an extrapolation of in- 
frared con tinuum or a power-law X-ray continuum, are included 
dSun & M alkan 1989; Laor 199G) (but also see Gaskell 2007). 

The standard thin disc is geomet rically thin, optically thick , 
and its emission is close to blackbody jShaku ra & Sunvaevlll973h . 
The gases in the disc rotate at Keplerian velocities, and they are 
in hydrodynamical equilibrium in the vertical direction. The struc- 
ture of a standar d thin disc can be derived analytically (see, e.g., 
iKato et aT]|l998h . In this work, we adopt the temperature distri- 
bution of an acc retion disc as a function of radius R given by 
iKato et al.l (l99^ . 
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Here M is the black hole mass, i?g = 2MG/c^, MEdd = 1.5 X 
10^*mgs~^,m = M/Mq, and a is the viscosity parameter. The 
spectrum of the accretion disc can be calculated with 
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where T(7?) is given by Eq. l[T), i?in = 3i?g, 7?out is the outer 
radius of disc, h is the Plank's constant, i is the inclination of axis 
of the disc with respect to the line of sight, and D is the distance 
from the observer to the black hole l lFrank et"al]|2002h . 

The optical/UV continuum of a standard thin disc is oc 
y^^^, w hich is incon sistent with fi, oc observed in many 
quasars. iGaskelll ( [2007,) suggested that the observed quasar spectra 
can be reproduced by accretion discs with a temperature gradient of 
T{R) oc 7?-° " instead of r(_R) oc R'"-"^^ as predicted by stan- 
dard thin disc model. Thus, we also adopt this modified disc model 
in our spectral calculations (this model is referred to as modified 
accretion disc model hereafter). The temperature of the disc is then 
given by 



T{R) = Ccor6.9 X 10 Q 
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where Ccor is a correction factor used to let the resulted disc lumi- 
nosity be the same as that for a standard thin disc with the same 
disc parameters. The spectra of the modified accretion discs can be 
calculated with Eqs. ^ and (O. 
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Figure 1. The spectra of standard thin discs/modified accretion discs for 
different black hole masses. The black and red lines correspond to stan- 
dard thin disc model and modified disc model, respectively. The solid lines 
correspond to mo = 0.1, while the dashed lines represent the cases with 
m = mo it <5m ((5m = 0.4mo). 



3 RESULTS 

The spectra of standard thin discs can be calculated with Eqs. (TJ 
and (|2}, provided the values of disc parameters (m, m and q) are 
supplied. In Fig. [T] we plot the spectra of the discs with different 
black hole masses and accretion rates. For comparison, we also plot 
the spectra of the discs calculated with Eqs. ([3} and ([2} based on 
the modified disc model described in last section. In the figure, we 
find that both standard thin disc model and modified disc model 
naturally show the optical/UV variability becomes larger with the 
increase of black hole mass, if the accretion rate changes at the 
same percentage. 

Assuming the optical variability to be caused by the change of 
accretion rates in AGN, we calculate the optical variability ampli- 
tude as a function of black hole mass with the disc models described 
in last section. In our calculations, the accretion rate varies in the 
range of m = mo ± (5m, where mo = 0.1 and (5m = 0.4rho for 
standard disc model ((5m = 0.5mo for the modified disc model), 
are adopted. In Fig. (2] we compare our calculations with the corre- 
lation betwee n variability amplitude at R-band and black hole mass 
discovered bv lWoldetal] ( [2007i) . We find that our calculations can 
reproduce their statistic result quite well(see Fig.|2](. 

We also calculate the optical variability amplitude as func- 
tions of luminosity by changing accretion rates m for given black 
hole mass m (see Fig. [3)- We find that the optical variability 
amplitude decreases with increasing luminosity for a given black 
hole mass, which is qualitatively consistent w ith the observed anti- 
correlation between these two quantities (e.g ., ICristiani et al .Il996l; 
IVanden Berk et alj2004IWilhite et al.ll2007l) . 

The dependence of the optical variability amplitude on the ob- 
served wavelength measured in the rest frame can be calculated 
based on either the standard thin disc model or the modified disc 
model for given disc parameters (see Fig. |4j. It is found that the 
variability amplitude decreases with increasing rest-frame wave- 
length is either for standard accretion disc model or the modified 
disc model. 
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Figure 2. The variability amplitude in R-band for accretion discs as func- 
tions of black hole mass. The black lines and red lines are for standard thin 
disc model and modified disc model, respectively. The dashed, solid, dotted 
and dashed-dotted line are for accretion rates mo=0.05, 0.1, 0.2 and 0.4 re- 
spectively. All the calculations are carried out by assuming accretion rates 
to vary with 5m ~ 0.4mo for standard thin disc model (<5m ~ 0.5rho 
for the modified thin disc model). The five points with error bars are the 
statistical results given in Fig. 5 of Wold et al. (2007) (see their paper for 
the details). The solid lines can fit the statistical results quite well. 



Figure 4. The variability amplitude of different accretion rate changes as a 
function of rest-frame wavelength. The arrangement of black line, red line, 
accretion rate and variation of accretion rate are the same with Fig.[T] 
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Figure 3. The variability amplitude as functions of luminosity for different 
black hole masses. The black lines and red lines are for standard thin disc 
model and the modified disc model, respectively, which are the same as in 
Fig.E 



4 CONCLUSIONS AND DISCUSSION 

The temperature of the disc decreases with increasing black hole 
mass, so the spectra of the discs surrounding massive black holes 
are softer than those for less massive black holes (see Fig. [Hi. It is 
found that the optical variability amplitude varies with black hole 
mass even if the changes of accretion rates rh are same. Our present 
calculations show that the observed correlation between the optical 
variability amplitude and black hole mass can be explained by ac- 
cretion disc models, if the variability is triggered by the variation of 
accretion rate. In all our calculations, a conventional value of vis- 



cosity parameter a = 0.1 is adopted. The final results depend on 
the value of a quite insensitively (see Eqs.[T]and|3j, and the con- 
clusions will not be altered even if a different value of a is adopted. 

In this work, the accretion rate m = 0.1 ± 0.04 (standard thin 
disc model) and m = 0. 1 ±0.05 (modified disc model) are adopted 
to calculate the disc spectra, and we find that the optical variability 
amplitude increases with black hole mass, which can reproduce the 
observed correlation quit e well (see Fig. [2l >. The mean accretion 
rate of the sample used in IWold et alj ( 120071) is ~ 0.0955, which is 
roughly consistent with the mean accretion rate m = 0.1 adopted 
in this work provided the radiative efficiency is ~ 0.1. 

The optical emission is dominantly from the accretion disc re- 
gion with several hundred Schwarzschild radii from the black hole 
(e.g.,,Gaskell 2007; Liu et al. 2008). The timescales of quasar vari- 
ability in lWold et alj ( |2007|) range from a few hours to years, which 
are much shorter than the viscous timescale of the disc, but is usu- 
ally longer than the light-trav el time of the distance to the black 
hole (see detailed discussion in lGaskelll2007 '). This may imply that 
the energy propagates at close to the speed of light, rather than on 
viscous timescales. Thus, the outer disc region is irradiated by the 
radiation of the inner disc region, and the temperature distribution 
of the disc is altered significantly as described by the modified ac- 
cretion disc model. Such a modified disc model can successfully 
reproduce the steep spectra in optical-UV bands l iGaskelli2007i) . 

Our calculations indicate that the optical variability amplitude 
declines with luminosity of the accretion disc if it black hole mass 
is fixed. This is qualitatively consistent with the statistical result 
given bv lWilhiteet"ai] ( l2007l) . Based on either the standard accre- 
tion disc model or the modified disc model, the spectral calcu- 
lation show that the optical variability amplitude decreases with 
increasing wavelength. Such a relation has already observed by 
Vanden Berk et al. (2004). 
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